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Abstract 

IT) 

In present article we consider the problems of concentrated point force 
, which is moving with constant velocity and oscillating with cyclic frequency 

in unbounded homogeneous anisotropic elastic two-dimensional medium. The 
properties of plane waves and their phase, slowness and ray or group velocity 
curves for 2D problem in moving coordinate system are described. By using 



the Fourier integral transform techniques and established the properties of the 
plane waves, the explicit representation of the elastodynamic Green's tensor 
is obtained for all types of source motion as a sum of the integrals over the 
finite interval. The dynamic components of the Green's tensor are extracted. 

The stationary phase method is applied to derive an asymptotic approx- 
imation of the far wave field. The simple formulae for Poynting energy flux 
^ , vectors for moving and fixed observers are presented too. It is noted that in 

the far zones the cylindrical waves are separated under kinematics and energy. 

It is shown that the motion bring some differences in the far field proper- 
ties. They are modification of the wave propagation zones and their number, 
fast and slow waves appearance under trans- and superseismic motion and so 
on. 



1 Introduction 

The theory of elasticity with moving source have a range of important applications. 
They are linked with the developing of more high-speed transport and with the 
necessity to evaluate the influence of elastic waves from moving objects on different 
constructions. The theoretical side of these problems is also of interest. The changes 
in the character of mechanical fields is essentially depend on the behavior of source 
motion and are caused by the type of differential equations ranging from elliptical 
to hyperbolic. 

From the steady-state problems with moving sources the problems with moving 
and oscillating sources are the most complicated. We use the following terminology. 



If the source is moving with constant velocity w and simultaneously oscillating 
with frequency w, we consider the problem B, while in the case w = 0, we have the 
classical harmonic problems with oscillating source and so refer them to the problem 



Nowadays considerable progress is obtained in the investigation of the problem 
B. The correspondence principles between the problems B and A [1], analogous 
known in fluid mechanic [2], are stated. The principles for unique solution selection 
are investigated [1, 3]. The energetic principles are formulated and the general 
theorems about energy transport are established [3], the some actual problems B 
for isotropic media are solved [4-6]. 

It is evident, that for anisotropic elastic media the particular problems B are not 
studied sufficiently In [7] the plane waves and fundamental solutions for problem 
B for tree-dimensional anisotropic media are studied, and in [8] the approaches for 
the problem B for anisotropic elastic and piezoelectric waveguides are proposed. 

In present paper we investigate the properties of plane waves and fundamental 
solutions in the problem B for anisotropic elastic plane. 

It is well-known, that for many classes of anisotropic materials it is possible to 
formulate the problems A and B in conditions of plane deformation. For existence 
of the plane problem several elastic modules must be equal to zero. Thus, for the 
plane deformation in the plane for elastic modules in two-index notation we 

require that C14 = C 34 = C 45 = Ci 6 = C 36 = C 56 = 0, for the plane deformation 
in the plane 0x±X2 Cu = C24 = C4Q = C15 = C25 = C56 = 0, and for the plane 
deformation in the plane OX2X3 C25 = C26 = C35 = C36 = C45 = C<±% = 0. 

The vectors of external force in the plane deformation problems must not have 
non-zero components in the perpendicular direction and must not depend on per- 
pendicular to plane coordinate. 

We assume that all above-mentioned conditions are implemented. To be more 
concrete, let us accept that the plane deformation is realized in plane OX1X2, and 
therefore the vector of mechanic displacement has non-zero components U\ and u^- 
u = {iii(x, t), «2( x ? t)}, x = {xi,x 2 }- For harmonic vibrations with frequency u> we 
shall find the established solution in the form 



In the problem with moving sources we consider two coordinate systems. Let 
C2} be fixed coordinate system with time r, and let {xi,x 2 } be coordinate 
system, which is moving relatively to fixed system with constant velocity w = 
{wi,w 2 }. We denote t the time in moving coordinate system. The two coordinate 
systems are connected with each other by the following relations 



For the problem B we propose that in moving coordinate system the harmonic 
behavior (1.1) with frequency u is existed. 

For the equations of theory of elasticity in moving coordinate system we use the 
following relations from (1.2) 



A. 



u = vexp (iut) . 




X! = £1 - WiT, X 2 = £2 - W 2 T, t = T. 



(1.2) 



V 5 = V X , d T = d t -w- V x , 



(1.3) 



where V x = {d/dx u d/dx 2 } = {d u d 2 }, V* = {d/d^d/dfr}. 

Then the stain vector S = {eu, e 22 , 2e 12 } {e^ are the components of strain tensor) 
may be expressed from the displacement vector as 

S = L(V X ) ■ u, (1.4) 

where 

L(V X ) = d 2 . (1.5) 
\d 2 d 1 J 

From Hook's low for anisotropic media for the plane deformation in the plane 
(xix 2 ) we have the relation between the stress vector T = {an, cr 22 , a\ 2 } (cr^- are the 
components of stress tensor) and the strain vector S 

T = C S, (1.6) 

where 

(Cn C12 Ci6 \ 
C 22 C 26 . (1.7) 
sym C 66 / 

We suppose here is not depended on coordinate, i.e. the medium material is 
homogeneous. 

By using (1.3) and the previous notations the equations for elastic medium in 
the case of plane deformation may be written in the form 

L*(V X ) -T + { = p(d t - w-V x ) 2 u, (1.8) 

where p is the density (p = const), f is the body force, (...)* is the conjugation 
operator. 

By using (1.4) — (1.6) the equation (1.8) may be rewritten 

L*(V X ) ■ C • L{V X ) ■ u + f = p(d t - w V x ) 2 u. (1.9) 

This equation together with the (1.1) is the equation for amplitude v. 
We study the properties of plane wave in problem B at first, to determine the 
fundamental solution. 



2 Plane waves and their characteristic curves 

We find the solution of equations (1.9) without body force (f = 0) in the form of 
plane waves 

u = Apexp [i(ut — ct ■ x)], (2.1) 

where A is the amplitude, p is the unit polarization vector (|p| = 1), a is the wave 
vector (a = an, |n| = 1, n is the unit wave normal vector). 

By substituting (2.1) into (1.9) yields the eigenvalue problem 

r(n)-p = pz/ 2 (n)p, (2.2) 



where 

v\n) = {c B p {n)+w n )\ (2.3) 

w n = w • n, Cp (n) = cj/q; is the phase velocity in the problem B, T(n) = L*(n) • C • 
L(n) is the acoustic Christoffel's tensor (matrix). 

As it follows from common properties of the problem A in 3D [9, 10], here in the 
problem A in 2D two plane waves with phase velocities z^(n) = c^-(n), c^-(n) > 0, 
j = 1,2, and with polarization vector (eigenvector) exist for arbitrary direction 
n,which can be selected orthonormalized. For the problem B, as it is obvious from 
(2.2), the situation is more complicated. Particularly, for the direct plane waves 
(Cpj > 0) in the problem B with fixed n and j G {1, 2} 

( Cj J(n) > w n ) A (<£.(n) > -w n ) =► c*(n) = 4.(11) - w n , (2.4) 

4(n) < -w n c*/(n) = (-l) fc 4(n) - «;„; fc = 0, 1, (2.5) 

Cpj(n) <«;„ =>- there are not direct plane waves. (2.6) 

Thus, in condition (2.4) for fixed n and j we have one plane wave, in condition 
(2.5) we have two (one - k — - fast, other - k — 1 — slow), and in condition (2.6) 
the direct plane waves are absent. The polarization vectors pj of plane waves in the 
problems A and B are identical, with the exception of the case (2.6), besides in the 
case (2.5) both fast and slow waves have the same polarization vector. Henceforth 
we denote the phase velocity in the cases (2.4), (2.5) by uniform way: = 

Cpj k \n) = (-l) fc c^(n) - w n , (k = 0) V (k = 0, 1), i.e. for the case (2.4) eg = c* (0) . 

If for Vn, Vj G {1, 2} the condition (2.4) is realized, then we shall call subseismic 
the motion behavior, and otherwise - trans- or superseismic. The difference between 
that notions we shall mark below. 

We introduce the phase velocity vector c^j^ and the inverse velocity vector 

C T = <f n = (-!)*<£■(*) - w nn; c&(n) = 4(n)n, (2.7) 

Lf (fc) = n/eff. (2.8) 

By analogy with problem A, we shall also call the vector refraction vector 

or slowness vector. 

The dispersion equation in the problem B for plane wave (2.1) have the following 
form 

D B (ot, u>) = det [r(a) - ptt 2 (a)E] = 0, (2.9) 
where E is the unit matrix 2x2, 

n{a) = uj + w • a. (2.10) 
The solution of equation (2.10) can be leaved in the form of the set of hypersur- 

fclCGS 

u = fc) (a) = «sf } (n), j = 1,2. (2.11) 



From the dispersion surface (2.11) the group velocity vector Cgj^ = Cgj k \n) 
may be found by the following way 



B(k) ® U j ( a ) B(k) . *x ® C pj ., ,., 

C <" = = °™ + ( E " nn ) ' ( 2 - 12 ) 



or, by using (2.7), 

cS W = (-l) fc c»-w. (2.13) 

It is known [10], that for plane waves in the problem A these important relations 
are correct (L A = n/c A ) 

c A -n = c A , cf.L A = l, n-dcf = 0, cf ■ dL A = 0. 

As it is shown in [7] for the case of tree-dimensional problem, the analogical 
formulae are valid. This proof is completely transported on the considered plane 
problem. Therefore the following relations are correct 

*T • n = c^; cff > ■ Lf « = 1; n • cfcff = 0; c£< fc > • dhf <*> = 0. (2.14) 

(In (2.14) and in the next the summation by repeating index is absent.) 

After analyzing obtained formulae (2.14), we may establish for the problem B 
the basic properties of phase velocity curves c^ k \n) = Cpj k \n)n, slowness curves 
(inverse velocity or reflection curves) L^ fe ^(n) = n/Cpj k \n) and group velocity 

curves (waves curves or ray curves) Cgj k \n). In general, all of these curves in the 
problem B with w ^ do not have the central symmetry n <-> (— n), and their 
crystallographical symmetry does not reserve. The phase velocity curves Cpj k ^ and 

waves curves are limited for all values w. But the slowness curves L^*^ can 
be both limited and unlimited depending on motion behavior. 

We shall call the source motion rate superseismic, if in R 2 the directions n, along 
which there are no slowness curves, exist, and all slowness curves are unlimited. If 
there is one limited slowness curve, but also unlimited slowness curves exist, then 
we shall call this motion behavior transseismic. 

We note that the slowness curves L^ k ^ may essentially differ from corresponding 
curves for the problem A, and their number may vary from 2 to 4. At the same 
time, the group velocity curves Cgj^ for the problem B may be obtained by simple 
transfer of curves for the problem A by vector (— w). Besides, with trans- and 
superseismic motion rate two parts and c^ 1 form one closed curve. For this 
behavior in case (2.5) for fixed n and j we have two plane waves (fast and slow) 
with inverse velocities, which belong to two different slowness curves. 

In addition to this from (2.2) and (2.14) important common properties follow 
[7]. Thus, the quasi-longitudinal and quasi-shear curves save their types for all w, 
because the polarization vector pj from (2.2) does not depend on the velocity w. 

The group velocity vector c g , by (2.14) in every point of slowness curve L, is 
orthogonal to the tangent to L (c s • dL — 0). Conversely, the wave normal vector n 
of plane wave with group velocity c s (n) is orthogonal to tangent in the corresponding 
point of wave velocity c g (n • dc g = 0). 



Everything considered above is illustrated by Figs. 1-11. These figures show 
slowness and group velocities curves for different kinds of quartz for plane deforma- 
tions in planes Ox\X 2 and OX2X3. Curves marked with "1" and "2" correspond to 
waves with phase velocities and respectively. The waves are numbered so 
that for the problem A < c^ 2 . Therefore in the most cases waves with subscript 
"1" will be quasi-shear and waves with subscript "2" will be quasi-longitudinal. 

The figures 1-6 show characteristic curves of plane waves for fused silica (SiO^), 
which is isotropic material. The elasticity modules and density of fused silica as well 
as other materials were taken from [9]. 

The figures 1 and 2 show curves of slowness and group velocities of fused silica 
without motion (w = 0) respectively. As it can be seen from Fig. 1 and 2 when 
w = these curves represent couples of concentric circumferences. 

When the velocity is not equal to zero pictures of characteristic curves essentially 
change. For subseismic motion when w\ = 0.6 • 10 3 m/s; w 2 = the curves slowness 
and group velocities are shown on Fig. 3, 4 respectively. Comparing Fig. 1 and Fig. 3 
we can notice that the curves of slowness change their shape and structure when w 
increases. Meantime, as it was mentioned above, group velocity curves (Fig. 2 and 
Fig. 4) simply transfer on vector — w relatively to the center of coordinate system. 
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Fig. 3. Slowness curves, Fig. 4. Group velocity curves, 

fused silica, w 1 = 2600, w 2 = 0. fused silica, w 1 = 2600, w 2 = 0. 
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Fig. 5. Slowness curves, 
fused silica, W\ = 7500, w 2 = 0. 



Fig. 6. Group velocity curves, 
fused silica, W\ = 7500, w 2 = 0. 



The changes of the shape and structure of slowness curves and ultimate ten- 
dencies of these changes can be seen from Fig. 5, 6, where characteristic curves for 
superseismic motion when wi = 7.5 • 10 3 m/s; w 2 = are shown. Here Fig. 5 shows 
slowness curves and Fig. 6 shows group velocity curves. Slowness curves for super- 
seismic motion (Fig. 5) are not limited. Each couple of these curves corresponds to 
fast and slow waves. Group velocity curves (Fig. 6) when wi — 7. ■ 10 3 m/s; w 2 = 
wholly move to the half-plane x < 0, this is caused by the absence of waves in front 
of the source for superseismic motion. The parts of fast and slow waves for group 
velocity curves are the pieces of the same curves separated by tangents to these 
curves which pass through the center of coordinate system (Fig. 6). 

The figures 5 and 6 illustrate relations of orthogonality c g ■ dL = and L • dc g = 
n • dc g = 0. Here slowness curves have directions n( ) and n( ) associated with one 
direction X( n ) on picture of group velocities, and vice versa. 

For anisotropic materials the pictures of characteristic curves can be much more 
complicate. The Figures 7-10 show curves for slowness and group velocities curves 
for a-quartz (Si0 2 ) for plane deformation in the plane x 2 x^. This material relates 
to trigonal crystal system of 32 class and possesses piezoelectric properties which we 

don't t a k p in^<^> rnnsirlpratinn 
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Fig. 7. Slowness curves, Fig. 8. Group velocity curves, 

a-quartz, w\ = w 2 = 0. a-quartz, w\ = w 2 = 0. 




Fig. 9. Slowness curves, Fig. 10. Group velocity curves, 

a-quartz, w\ = 8000, W2 = 0. a-quartz, w\ = 8000, W2 = 0. 

Slowness curves for a-quartz on Fig. 7 and 9 have the points of inflection and 
the sections of convexity and concavity. Therefore group velocity curves have typ- 
ical acute edges. Availability of such edges defines areas with different numbers of 
spreading waves even for problem A [11]. This property is also secured for problem 
B with additional sophistication of the wave picture and distraction of the curve 
central symmetry. 



3 Fundamental solutions 

The fundamental solution for the problem B is the solution of equation (1.9) with 
the point external source f 

f = /5(x)l exp (iuot). 

For obtaining unique solution we use the principle of limiting absorption, ac- 
cording to which in the case of oscillation of point source by the law exp (iut) the 
amplitude v(x) in equation (1.1) should be defined as a limit at e — > +0 

v = lim v e . (3.1) 

e^+0 

The function v e satisfies the following e-problem 

-L*(V*) • C • L(V X ) ■ v £ + p{ioo £ - w • V x ) 2 v e = /<5(x)l, (3.2) 

where uj £ — uj — is, e > 0, e <C 1. 

By applying Fourier integral transforms along xi, X2, we obtain the fundamental 
solution v e for equation (3.2) in integral form 

v e (jc) = J- ff K _1 (a) -le-^daidaa, (3.3) 
4vr 2 J J R 2 

where 

K(a) = r(a) - p^(a)E, (3.4) 
£l £ (ct) = uj £ + w • a. (3.5) 



Introducing the polar coordinate system (a, 6) 

ol = an, n = {cos 9; sin#}, 

it is possible to decompose the matrix K _1 (o:) by accompanying matrices P 
which are composed of eigenvectors pi, p 2 for problem (2.2) 



P, = P J (n)=p,(n)p*(n), n = n(9), 
with regard to (3.6) — (3.8) we obtain from (3.3) 

2 ,2tt 



/ 



.7=1 J ° 



+oo 



Ij£ ~ J aVf (n) - Q2(cm) 6 da ' z ~ n x - 
We decompose integrand function from (3.10) by partial fractions 

1 / at a - 



aV 2 (n) - f^( an ) 2iv e Vj I a - at a - a je I ' 



Us _ U £ 



J 7/. — in J 



Vj - W n J Vj + W n 

In the integral (3.9) from 9 we realize the following transformations 

2?r tt/2+8 tt/2+§ 

j {...){&) d9= j (...)(6)d0+ j (...)(6 + n)d0, 
-tt/2+9 -71-/2+0 

where 9 is the angle in the polar coordinate system for physical plane x 

x = ry, |y| = 1, y = {cos 6, sin9}. 

By taking into account (3.11) — (3.13) the formulae (3.9), (3.10) can be 
formed in the following way 

2 

v £ (x) = :-^— E / ^_lj J£ (9)d9, 

J -tt/2+9 

/ r+oo g—iza r+oo ^iza \ 

Jj £ (Q) — at ( / — da+ T da - 

3£ \J a- aj £ J a + aj £ j 



r+oo Q—iza r+oo ^iza 

—aj £ I / — da + — da 

V Jo ol- a j£ J a + a j£ 



The calculation of integrals from (3.16) depends on real and imaginary parts of 
t , which in their turn depend on moving rate 

Rea^ £ > 0, Imat < 0, Vj — w n > 0, 



Re at < 0, Imat > 0, Uj -w n <0, (3.17) 

Re<x~ < 0, Ima^ > 0, Vj + w n > 0, 

Re > 0, Im aj £ < 0, Vj + w n < 0. 
Besides, it is essentially, that V6> e [-7r/2 + 7r/2 + 6>] 

,2 = n-x = rcos(#-#) > 0. (3.18) 

We shall introduce the contours Cp = [0, R] U U [±ii?, 0], where and 
are the quoters of circle with radius R and center in the origin of coordinate system, 
lying respectively in the first and in the fourth quadrant of complex plane a, R >> 1. 

Calculating the integrals by complex integration method, using contours Cp , 
relations (3.17), (3.18) and Gordan's lemma and turned R to infinity, we obtain 
(z > 0) 

-+oo — iza 



± r + °° [a — in) 

da = -l-KiQ-^H^V; -w n ]+i / -—^5 — A:e~ zv dr], (3.19) 

Jo (Q!jJ +V 

,+oo e izg ,+oo ( ± _ ^) 

/ T ± da = i ' e 77 cfy, 

where if is the Heaviside function. 

On substitution (3.19) into (3.16) and received formula into (3.15) and imple- 
menting passage to the limit when e — > +0, we obtain the final representation of 
fundamental solutions for the problems A and B with arbitrary moving rate 

v = w d + v , (3.20) 



ir/2+6 



/ e ^ (3 ' 22) 



" ( W> 



9 tt/2+0 



' = ~-tt/2+6> 



= / \ a >T^f- - a '^0? 1 e ~" (3 24) 

z = n • x = r cos(# — 9), af = uo /(±Vj — w n ), (3.25) 

where symbol Yl'k i n (3-21) denotes the presence or the absence of summation by k 
according to moving rate. 

In addition formula (3.12) can be transformed by using the following integrals 

r+oo -zt) 

/ 7J— ~ Vdrj = - [ci (C*) cos(C^) + si (C*) sin(C^)] , (3.26) 
Jo C 2 + V 



/ a 2 , 2 dr l = 7 t ci (C 2 ) sin (C^) - si ((z) cos(C^)] , 
Jo C + T C 

where si (x) and ci (x) are respectively integral sine and cosine. 

For the problem A with w = the fundamental solutions representations in form 
(3.21) — (3.26) can be significantly simplified 



... n/2+e 



" = £E / ^ 

3 -w/2+e 

2 tt/2+9 

v o = -TT^Z^ / C1 cos + si -z sin -z dg. 

J -n/2+§ 

Thus, we have obtained the fundamental solutions for problem A and B in 
the form of integrals by finite interval. It should be noted that the representation 
(3.20) — (3.26) for the problem B are suitable for arbitrary moving rate. 



4 Far field asymptotics 

The term from (3.21), (3.22) or (3.27) define dynamic effect of displacement field. 
We shall determine the asymptotics of in far field our >> 1 for general case of the 
problem B. For this it is evidently required to find the asymptotics of oscillating 
integral (3.22). 

From classical method of stationary phase the contribution of individual station- 
ary point 9 S for integral 

J F(9)e iq( - e)ujr d9 (4.1) 
at ujr » 1 can be given by the expression 

, 27r , Ffflje^^+i 8 ^ 9 "^ dO, (4.2) 
ur\q"{9 s )\ 



where 9 S is the saddle (stationary) point, or the root of equation 

q'(9 s ) = 0. (4.3) 
In our problem for fixed j and k we have 

We shall obtain the formula for saddle point 9 S . From (4.3), (4.4) we have 

(4.5) 



dq 



09 89 " 

We can represent the components of vector dLf^/d9 in the form 



and since 



it follows 



09 ^ dn m 09 ' 



dn rn dn rn \^) } c Bf) K ° im ^ dn rn h 
dL f _ 1 r F jB {k)( d$ k) dn 



Then by using (2.12) we find 



B(*) _ <9cg fc) 9n 
9j 89 dn 89' y } 



since • (dn/09) = 0, n • (dn/90) = 0. 

From (4.7) the formula (4.6) can be rewritten in the form 



FIT B ( fe ) 1 Q 

1 -[E-LfWri-t- (4-8) 



00 B(fc) L 3 ^ 93 QQ 

U P3 



On substitution of (4.8) into (4.5), (4.3), we obtain 



y* • [E — Lf»(«£«)1 -|5 = 0. (4.9) 



Besides, using one of the relations (2.14) (L? ( k) ■ c B g ^ k) = 1) we can write following 
sequence of equalities 

= n.y-(Lf fe ).cff))n.y = n*.E.y-(n.cff)^.^- 



93 I J J \ 93 I B(k) 

C P3 



= n* ■ E ■ y - n* • (cf/ fc) (Lf (fc) )*) • y = y* ■ [E - Lf (fc) (c* (fe) )*] ' n, 



y*.[E-Lf (fc) (cf/ fc) )*]-n = 0. (4.10) 



93 ^ 3 / J •> J L j \ gj 

i.e. 

v * . [F, — ~\_ . 

3 \ 93 

From (4.9), (4.10) we obtain, that the vector y* • [E-Lf (fc) (cJ (fc) )*] is orthogonal 
for two one-to-one orthogonal vectors n and dn/d6, and therefore, it is equal to zero 

y *.[E-Lf fc) (c* (fc) )*] = 0. (4.11) 
By using inequality (3.18), which can be rewritten in the form 

y.Lf fc) (0)>O, (4.12) 
we can get the following shape for the equation (4.11) 

4 {k) (6)\yLf {k \6)\=y. (4.13) 

The relations (4.13), (4.12) are the suitable formulae for definitions of stationary 
points 9j^l and wave normal vector = n(9j2). 

Geometrically, from (4.12), (4.13), just as in problem A, for stationary points the 
group velocity vector Cgj^ = c^ fc ^(n^), which is perpendicular to slowness curve 
hf^ k \ is directed along vector y (x). Therefore, for fixed direction x on physical 

(k) 

plane the stationary value will appear such values n, in which the external 
normal to slowness curve coincides with direction x. 

Now, we shall dwell upon calculation the expression q" in the stationary point 

8 2 L B(k) 

q " = ~ y • ~w^- (414) 

Let r be tangent to slowness curve L^*^ unit vector and s be natural parameter 
along this curve. We note, that y is the unit normal vector to ~L B ^ k \ and thus, 
ydLf (k) /d9 = 0. 



3 

Because 



dy dL? k) _d_( dLf w \ ^ PL?™ 



39 89 89 r 89 ) J 89 2 ' 
then we can transform formulae (4.14) in the form 

By using Frene formulae, we can write 

8Lf^ = 8Lf^8 1= 8s_ 8y = 8yd l= B{k) 8s_ 
89 8s 89 T 89' 89 8s 89 pj T 89' 1 ' ' 

where is the curvature of slowness curve h B ^ k \ 



From (4.15), (4.16) we have 



„__ k B(k) w\ _ jB(k) ""J ( . 7) 

V - K pj \ xa I - K P.i xa xa ■ V^- 1 ') 



'ds\ 2 _ k s {k) d_L^_ 9Lf {k) 
'"> \06) ~~ pj 06 ' 06 

By using (4.8), it is easy to find, that 

06 06 (c*} k Y PJ 93 d6 )V 1 ' 

By decomposing group velocity vector c B j k ^ in system of orthonormalized vector 
n, On/ 06, we have 

B(k) _ B(k) B(k)dn 
^gj u gjn ^ u gjd qq > 



with c $> = <£<*>, c $> = ■ (On/06). 
Therefore, 



I B(k)\2 _ (B(k)\2 , ( B(k) 9n\ 
l%i \ u p.i > " r Sm xn / ' 



' \ ~gj qq J 
and as the result we obtain from (4.17), (4.18) the simple expression 

„ k B(k) ^ gj ) (419) 

By substituting of (4.19) into (4.2), (4.1), (3.27) yields the following formulae 
for far field asymptotic 

v- = £(£)'£*$, ( 4 - 2 °) 

j=l k m=l 

V (fc) » (-l) fc+1 • 1 Sj™ e - <(hff . L ?W. y+f ff ?W) (4 

Jm o /o — i B ( k h a f\l B ( k )\ 

zy zirr p\Cgj m \c pjm y \K pjm \ 

where is the number of stationary point for separated slowness curve, ~L B ^ = 

jB(k) ( (k), B(k) = B(k) ( (k)\ B(k) B(k) ( (fc)v 4 . . „A/ n (fc)x r.S(fc) . . 

j V jm)i gjm gj \ jm)i pjm pj \ jm)i ^pjm ^pj\ ll jm) ' ^pjm 

k B(k) (n {k) ) a B{k) = Hienk B(k) 

pj \ jm)i jm o pjm ' 

Because of (— l) fc c^^/c^ m = u;/£)(a^), and in the stationary points from 
(4.12), (4.13) L^j^ • y = |c^^] _1 , then we can rewrite formula (4.21) in the form 



v<*> « (-!)*+! ^ (n ^ ) - 1 e " l te +fCTr /. ( 1.22) 



2v^pfi(af2)|c^ ) |^ F 



As we can see from (4.21), (4.22), in far zone the wave field is separated in 
individual cylindrical waves j = l,2;/c = 0,lorA; = 0;m = l,..., Nj k \ 



We note that the wave fields in the far zone have peculiarities in the neighborhood 
of directions n^(x) for which k^ m ~ for near located and for multiple points. For 
such directions as in the problem A [13] another shapes of asymptotic decomposition 
than (4.20)-(4.22) are required. 

As it follows from (4.12), (4,13), the number of waves vj^ in far field are de- 
termined by the number of stationary value n^(x). As we marked above, in the 

stationary value n^(x) the group velocity vector c^^(n^) is directed along the 
vector x. Therefore, for fixed direction x the number of waves is easy evaluated 
as the number of intersection between the ray Ox and the group velocity curves. 
Depending on moving rate, the number of waves may be essentially changed, and 
besides in the case of trans- and superseismic motion the zones of a fast and slow 
waves propagation, limited by Mach's cones, exist. 



The Figures 11, 12 show typical pictures of group velocity curves in planes 
which correspond to the work cartesian coordinate system. Separately the selected 
directions x on the Figs. 11, 12 are supplied with arrows and marked with numbers 
in circles which denote the number of waves for these directions, i.e. the number of 
intersection points of the curves cfj and Cgj k ^ with directions x. On the figures the 
points of intersection are bold-face points. Figs. 11, 12 illustrate, that the number 
of waves in different zones takes any values from to 8. 

In pole coordinate system (r, 8) the components of group velocity vector in 
stationary points are given by: c^ mr = |c^ m |; c B g - m Q = 0. Hence, cylindrical waves 
(4.21), (4.22) in far field satisfy all general conditions for cylindrical waves from [3]. 
Therefore, for waves (4.21) or (4.22) the general energetic relations are correct [3], 
and the group velocity vector c^ m is equal to energy transport velocity vector (ray 
velocity vector) for fixed and moving observers. 

We use general formulae [3] for the cylindrical waves average energies (E^ x ) and 
(E x ) in moving coordinate system for fixed and moving observers respectively 




Fig. 11. Group velocity curves, 
a-quartz, w\ = 5500, W2 = 



Fig. 12. Group velocity curves, 
a-quartz, w\ = 8000, u>2 = 0. 




(4.23) 




(4.24) 



and formulae for average energy flux vector (Poynting's vector) 



(jfW) = (4.25) 

where (...) = ± / Q T (...) di, T = 2vr/cu. 

By substituting of (4.22) into (4.23) — (4.25) yields the following expressions for 
energy flux of individual cylindrical waves in far field 



/ t£..-U-)\ ^ 2 / 2 | P i( n jm) • jj , , ; 

W . mr ; _ 



B(fe) 

/ T^(fe)\ _ U I T&{k)\ _ / 1 \k C Pjm I j£x{k)s /, 9? n 

\ J jrnrl — (i), W]tnr / — V ± / „A W 3 'mr /■ V^ ,Z ' / 

For superseismic moving rate fi(a] m ) < for k — 1, and therefore (J^ r ) < 0, 
i.e. the slow waves transfer the negative energy, measured by moving observer. But 
this property is common for slow waves in problem B with superseismic moving 
sources. 

The Figures 13-22 show the average energy flux curves (Jr X ^} for fused silica 
and a-quartz (Si0 2 ). The normalizing factor ( for the values of average energy flux 
is equal to uo 2 f 2 / '(IQtt pr) . 

The energy flux curves have more complicated forms, then slowness and group 
velocity curves. Really, according to (4.26) the values of energy flux depend on the 
values of group velocity, curvature of slowness curve and form the scalar product 
between polarization vector = pj (n^) and the unit vector of source direction 1 

(by (3.8) |Pj(n^) • 1| = |pj(n^)|). For some angle 9 the polarization vectors can be 
orthogonal to unit vector of source direction 1, because quasi-longitudinal and quasi- 
shear polarization vectors are perpendicular, and if the angle 6 change by 2n, the 
polarization vectors also rotate by 2ir. In these cases the energy flux vector is equal 
to zero. On the other hand the energy flux tends to infinity, when then curvature of 
slowness curve tends to zero. This exists, when slowness curves stretch on infinity 
and in the points of inflection. As a result, taking into account the possibility of 
existence of zones with different number of spread waves, the energy flux curves are 
more complicated and complex for analysis. 

The Figures 13-18 show the results of energy flux calculations for isotropic 
material fused silica. The values of source motion are equal to corresponding values 
for Figs. 1-6 (for Figs. 13, 14 w\ = w 2 = 0, for Figs. 15, 16 w\ = 2.6 • 10 3 m/s, for 
Figs. 17, 18 w\ = 7.5 • 10 3 m/s). The consistent pair of the figures correspond to 
directions of unit vector of source along motion direction (l± = 1, l 2 = for Figs. 13, 
15, 17), and perpendicular to motion direction (/i = 0, l 2 — 1 for Figs. 14, 16, 18). 
As it was used above, the curves marked subscript "1" correspond to the results for 
quasi-shear waves and the curves marked subscript "2" correspond to the results for 
quasi-longitudinal waves. 



Fig. 13. Energy flux curves, 
fused silica, W\ = w 2 = 0, li = 1. 



<jf >M0" 8 
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Fig. 14. Energy flux curves, 
fused silica, w± — w 2 — 0, l 2 — 1. 



<jf >M0" 8 
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Fig. 15. Energy flux curves, 
fused silica, W\ = 2600, w 2 — 0, l± — 1. 



<js x >,mo- 




Fig. 16. Energy flux curves, 
fused silica, W\ = 2600, w 2 = 0, / 2 = 1. 
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Fig. 17. Energy flux curves, Fig. 18. Energy flux curves, 

fused silica, wi = 7500, w 2 — 0, h — 1. fused silica, w\ = 7500, w 2 = 0, l 2 — 1. 

From the Figs. 13, 14 it is obvious, that in the problem A for isotropic material 
fused silica there is the maximum of energy flux for quasi-share wave in direction, 
which is perpendicular to direction of unit vector of source. All curves of energy 
source in the problem A are symmetrical to the origin of coordinate system. For the 
problem B with subseismic motion (Figs. 15, 16) the curves of energy source change 
significantly and become not symmetrical. For superseismic motion (Figs. 17, 18) 
the energy flux curves are not limited, because not far from boundaries of wave 



propagation (Mach's cones) the curvature of slowness curves tends to 0. As for the 
case of slowness curves each couple of energy flux curves with superseismic motion 
corresponds to fast and slow waves. 

The figures 19-22 illustrate the pictures of energy flux behavior for anisotropic 
material (a-quartz). For the Figs. 19, 20 we have the case of problem A with w = 
(analogically Figs. 7, 8), and for the Figs. 21, 22 we have the case of problem B with 
superseismic motion (analogically Figs. 9, 10). In this cases we obtain exceptionally 
complicated behavior of energy flux vectors. (For obviousness the subscript "1", 
which corresponds to quasi-shear waves is absent in Figs. 19-22.) 





Fig. 21. Energy flux curves, Fig. 22. Energy flux curves, 

a-quartz, w\ = 8000, w 2 — 0, l± — 1. a-quartz, w\ = 8000, w 2 — 0, l 2 — 1. 

In the conclusion we announce similar paper [14], where additional pictures of 
phase velocity are presented, but the results of energy flux calculation are absent. 
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